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Abstract 

The controlled growth of organic layers on surfaces is still waiting for an in-situ reliable technique that would allow their 
quality to be monitored and improved. Here we show that the growth of a perylene monolayer deposited on Ag(llO) at 
room temperature can be tracked with low energy atoms in a regime where the energy perpendicular to the layer is less 
than 0.1 eV and below the organic film damage threshold. The image processing required for this atom triangulation 
technique is described in detail. 


1. Introduction 

Thin films of polycyclic aromatic hydrocarbons such as 
perylene (C 20 Hi 2 ) are investigated intensely for their in¬ 
tegration in a broad range of applications, for instance 
in organic light emitting diodes (OLED). A large num¬ 
ber of studies over the last decade have demonstrated that 
molecular organization is a crucial factor for the physi¬ 
cal properties of organic layers. Thin film assembly can 
be investigated locally in great detail by STM P and/or 
AFM [H but not during the growth of the layer. Other 
techniques such as reflection high energy electron diffrac¬ 
tion (RHEED) widely used in molecular beam epitaxy, 
are used to probe growth in-situ but the high beam en¬ 
ergy destroys the fragile organic layers. In fact, there is 
a clear lack of techniques able to monitor in-situ and in 
real time the structure of the molecular network and its 
quality. Using the same grazing incidence geometry as 
RHEED, i.e. compatible with in situ operation, two tech¬ 
niques sensitive only to the topmost layer have emerged 
recently. These are grazing incidence fast atom diffrac¬ 
tion (GIFAD or FAEU Ij, 1^ and fast ion or atom beam 
triangulation (IBT) [j, y, @, 0, [1| • The first one is very 
demanding in terms of long range order while the second 
one, developed for inorganic crystals operates at compar¬ 
atively high interaction energy with the surface. In this 
paper, a new regime of extremely weak interaction with 
the surface is briefly presented, and applied to the investi¬ 
gation of perylene assembly on Ag(llO). A simple image 
processing method providing a more efficient triangulation 
technique is described in detail. The convergence between 
GIFAD and IBT is discussed at the end. 
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1.1. IBT and GIFAD 

IBT is a technique allowing the identification of the low 
index directions of the topmost layer of a crystalline sur¬ 
face. The working principle of this method is that an 
atom or an ion impacting on a surface at a grazing in¬ 
cidence angle will experience a specific physical condition 
known as axial channeling if it is aligned with a low index 
direction. Several physical properties can be exploited de¬ 
pending on the nature and energy Eq of the projectile. 
Protons with £'o=25 keV at an incidence angle ^^^=1.6 ° 
with respect to the Gu(OOl) surface have a normal energy 
E± = Eosin{0in)‘^ ^ 20eU. Monitoring the electric cur¬ 
rent flowing to the target during the azimuthal scan is a 
simple technique giving an indirect measure of the emis¬ 
sion of secondary electrons In this energy range and 
along the surface channeling direction, the projectiles un¬ 
dergo multiple bounces @ on the atomic rows (the walls 
of the channel) and have more chance to trigger electron 
emission. This interpretation was confirmed using an elec¬ 
tron detector in front of the surface, which provided a 
higher sensitivity allowing the beam intensity to be re¬ 
duced 0. At such comparatively high projectile normal 
energy, the initial projectile charge state has little con¬ 
sequence and similar contrast was observed with neutral 
hydrogen atoms 0. Only recently has the same trian¬ 
gulation strategy been applied to lower energy atoms of 
2 keV with a perpendicular energy E_\_ between 0.6 and 
1 eV 0, g. In this case, there is no electron emission 
and the azimuthal scan is performed by monitoring the 
variation of the intensity of the scattering profile in the 
specular region. The azimuthal contrast is generated by 
the fact that along the channeling directions, successive 
momentum transfers are cumulative instead of being ran¬ 
domly distributed. As a result the scattered beam has a 
lower, broader peak intensity with dips generated at each 
low index direction. In this range where the perpendicular 
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energy is only few meV, the distance between the projec¬ 
tile to the surface layer is around 1-3 A from the top layer. 
This is true only if the deposited layer is already well or¬ 
ganized or dense enough, so that mutual shadowing can 
occur via successive deflections which prevent the projec¬ 
tile from approaching too close to the molecules. On the 
other hand, if an isolated molecule or a terrace edge is en¬ 
countered on the bare metal surface, then the interaction 
energy is probably large enough to destroy the molecule 
or to penetrate the terrace. Here, we explore the same 
technique at even lower normal energies, with 300 eV he¬ 
lium atoms having an in the 50-100 meV range where 
the minimum distance to the surface is above 3-4 A so 
that even an isolated molecule will not experience a close 
collision, provided that it is lying flat on the surface. The 
helium atom will fly over the molecule, or the terrace edge, 
without penetration and with limited momentum transfer. 
Initially, this energy interval was chosen because it corre¬ 
sponds to the conditions where diffraction was observed 
in GIFAD. If the transverse coherence of the atomic 
beam is larger than the lattice parameter d, diffraction 
can be observed as Bragg peaks separated by a recipro¬ 
cal lattice vector G = 27r/d in the scattering profile. The 
transverse coherence is a simple statistical property of the 
beam expressed as = 27r/A/c with A/^ = koAip, where 
ko is the momentum of the atoms, and Aip is the angu¬ 
lar divergence of the beam. The additional information 
provided by the lattice parameter is extremely valuable to 
complement the triangulation. Indeed, triangulation pro¬ 
vides the directions of alignment but does not specify ei¬ 
ther the unit cell size or the number of molecules in the cell 
whereas even weak diffraction is sufficient to obtain these 
parameters with a quite high precision. In the case of a 
surface exhibiting sharp distinct diffraction spots, the pro¬ 
file of the surface electronic density can be retrieved with 
exceptional accuracy Eini- However, it should be noted 
that this high resolution regime has not been achieved so 
far for molecular layers. Clear diffraction lines have re¬ 
cently been reported for comparatively small, non-planar 
alanine molecules chemically adsorbed on a Cu(IIO) sur¬ 
face ii. The present paper reports the experimental 
investigation of perylene growth on Ag(IIO). This large 
planar molecule shows significantly lower interaction with 
the surface leading to a partial flexibility in the assembly 
of such an organic monolayer Q. Only weak diffraction 
is observed with a contrast weakening over time. This is 
probably due to the metastability of the layer. To comple¬ 
ment and prepare the diffraction measurement, we have 
developed a new simple and robust procedure for triangu¬ 
lation with low energy atoms. 

2. Experimental setup 

The deposition of a perylene monolayer on the Ag(IIO) 
surface was conducted in an ultra-high vacuum (UHV) 
system equipped with standard facilities for sample prepa¬ 
ration and analysis as well as a separate load-lock chamber 



Figure 1: Scheme of the experimental setup. The surface is placed 
on a Omicron plate depicted in gray in the center and exposed to 
the 300 eV He° beam at grazing incidence. The scattered beams are 
detected almost a meter downstream on a position sensitive detector 
imaging the outgoing angles. The primary beam and the detector 
remain fixed during the experiment. 


compatible with a high performance portable UHV suit¬ 
case. The Ag(IIO) surface was cleaned by several cycles 
of Ar ion sputtering followed by annealing. The perylene 
molecules were evaporated on the Ag(IIO) surface held at 
room temperature, and the molecular structure was sub¬ 
sequently checked with a microchannel plate-low energy 
electron diffraction (MCP-LEED) system operated at an 
electron energy of 18 eV, which is surface sensitive and 
limits the damage to the organic layer. The sample was 
then transferred by means of the portable UHV suitcase 
into the independent GIEAD chamber without breaking 
the UHV conditions. The surface is mounted perpendic¬ 
ular to the rotation axis of a standard X,Y,Z,(/) sample 
manipulator itself mounted on a custom designed tilt sys¬ 
tem controlling the angle of incidence 0 (see Eig. [T]). The 
atomic beam of 300 eV neutral helium passes into the vac¬ 
uum chamber through a 100 jim circular hole. Prior to 
inserting the sample, the beam position is recorded onto 
the position sensitive MGP detector in front of a phosphor 
screen and filmed by a GGD camera. Then the surface 
is progressively inserted in the beam and the angle of in¬ 
cidence is tuned to 0.7°. Most often, the surface is only 
partially inserted into the beam leaving a small part of the 
neutral helium atoms flying over the surface without in¬ 
teraction to produce a small reference spot on the detector 
as can be seen at the bottom of the scattering patterns in 
Fig. H 


3. Results 

A diffraction signal was observed after the sample was in¬ 
troduced into the beam and vanished few hours later under 
exactly the same conditions. Therefore, a systematic az¬ 
imuthal scan was performed. In principle, such a scan re¬ 
quires only a continuous or stepwise rotation synchronized 
with a movie or with a sequence of images. This was not 
possible because the surface plane defined by its normal 
vector was found to be misaligned with the rotation axis 
by approximately 1°, as can be seen in the azimuthal scans 
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Fig. [6]below. An automatic compensation could have been 
programmed as explained in [l2|. In the present case, ap¬ 
proximate manual corrections of the angle of incidence 0 
were needed every few images to maintain the scattering 
profile in the same region on the detector. Each image 
of the scattering profile is composed of fifty exposures of 
2 s accumulated directly inside the CCD camera before 
transfer to the PC. Consequently, in order to perform the 
azimuthal scan in a reasonable time, comparatively large 
angular 0 steps of two degrees were chosen. In such a 
long exposure, the camera produces a significant noise, so 
a dark image was also recorded under exactly the same 
exposure conditions but with the beam valve closed. This 
dark image is automatically subtracted from all images. 


3.1. Data analysis 



Figure 2: Raw images of the scattering patterns recorded every two 
degrees azimuth for a perylene monolayer on Ag(llO) substrate. The 
relative orientation is as sketched in Fig. ^ i.e. the surface plane is 
horizontal while the outgoing scattering angle 6out increases along 
the vertical direction. The color scale is a fixed black and white 
banded palette highlighting a quasi-elliptic white contour. The cor¬ 
responding labels (a) to (f) are reported on the azimuthal scan in 
Fig.m 

Six successive images of a region of interest in the de¬ 
tector are displayed in Fig. At first glance these im¬ 
ages look rather similar having a smooth scattering profile 
without any evidence of Bragg peaks. A closer inspection 
reveals that the pattern experiences profound variations. 
The white quasi-elliptic contour line undergoes an anti¬ 
clockwise rotation between the images in Fig. [Ifa) and 
Fig. [2jc) whereas the evolution between the images in 
Fig. [2fd) and Fig. [2]^f) is clockwise. Between Fig. [2]^c) 
and Fig. [2jd), the variation is more drastic with a sudden 
mirror change of the orientation. This behavior is charac¬ 
teristic of surface channeling of ions or atoms along a low 
index direction. Since the energy is more than 100 times 
less than in Q we try to clarify the specific dynamics oc¬ 
curring as the projectile direction approaches a low index 
direction as a function energy. 

At high normal energy E± ^10 eV the projectiles pene¬ 
trate deeply the surface electronic density. The curves of 


iso-electronic density no longer form a 2D manifold over 
the atomic layer but split into separated spheroids centered 
around the atoms so that some projectile trajectories may 
even penetrate below the surface. However, most of the 
projectiles simply bounce in-between the rows of surface 
atoms [l^ and undergo multiple collisions from one side 
of the channel to the other in a genuine guiding that gen¬ 
erates specific rainbows j^ . 

At E_\_ < 100 meV used in this work, the evolution as 
the beam direction approaches a low index direction has 
been measured in detail in the diffraction regime dllli 
and is well described by quantum mechanics Il6l or semi- 
classical models for weak surface corrugation [l7|. In this 
regime the net deflection of the beam is very limited and is 
localized at scattering ang les significantly larger than the 
specular angle [IJ, [l^, ll6| . In the present experiment the 
later effect is likely to be responsible for the clear rotation 
of the top of the contour ellipses in Fig. [2j As noted in 
[l5| this indicates that the twist angle of the ellipse is a 
differential signal which has a zero twist when the beam 
is perfectly aligned with a low index or along a random 
directions. To identify the low index directions, we con¬ 
centrate on robust features of the scattering pattern such 
as the horizontal or vertical intensity profiles. The associ¬ 
ated widths of the observed distribution seem an obvious 
choice, however, we have to take into account that the an¬ 
gle of incidence Oin could not be kept constant due to the 
tilt angle with the manipulator axis. We decided to trans¬ 
form our 2D images into a coordinate system which is less 
sensitive to the variation of the angle of incidence. 


3.1.1. Polar coordinates 

The Cartesian coordinates on the raw images in Fig. [2] 
correspond to scattering angles parallel and perpendicu¬ 
lar to the surface. For scattering problems in general, the 
magnitude of the deflection and the associated polar pro¬ 
file are the most relevant. Examples include (i) scatter¬ 
ing rainbows observed at at fixed polar scattering angle a 
(see Fig. 13) over a wide range of incidence angles and (ii) 
supernumerary rainbows [181. Polar coordinates are still 
valid in the quantum scattering regime where diffraction 
charts show clear nodal structures well-aligned with the 
polar angle. In fact these nodal structures transform into 
attenuating supernumerary rainbows when the diffraction 
is no longer visible. The use of these coordinates should 
provide image properties that are much less sensitive to 
the actual angle of incidence. 


3.1.2. Polar-like 2D transformation 

A pragmatic choice is to consider the position of the direct 
beam as a reference in the raw image as well as for the 
transformation. Taking (0,0) for the beam coordinates, 
each pixel of coordinates {kx^ky) is part of a circle having 
a center {kx^kyc) given by kxc = 0, kyc = ky/2 + kl/ky, 
a radius keff = kyc and a = sin~^{kx/keff). This polar- 
like transform (kx^ky) {a,2keff) is defined everywhere 


3 

















above the direct beam {ky >0). In order to provide a mea¬ 
sure close to the original one, the polar angles are reported 
between -90 ° and -1-90 ° only, corresponding to the upper 
parts of the circles in Fig. [3l Note that this transforma¬ 
tion does not affect at all the vertical line above the direct 
beam, i.e. the scattering profile in the incidence plane. 
Also, it does not seriously impact the general shape of the 
scattering patterns as seen in Fig. [3l 



-90® 0® +90® 


Figure 3: Sketch of the polar-like transformation associating (a) a 
coordinate (kx, ky) to (b) new coordinates (a, 2keff) defined by: the 
radius of the circle k^ff, the point {kx,ky), its left/right symmetry 
{—kx^ky) , and the direct beam. This forces the center to lie on the 
vertical axis above the primary beam, i.e. in the plane of incidence. 
Finally, a is the polar angle along this circle. In this transformation, 
the incidence plane (vertical line) is strictly invariant and most of 
the image is hardly affected. 


3.2. Azimuthal scans 

Fig. [5] displays a detailed azimuthal scan of the measured 
Lorentzian width Wa evidencing the presence of low in¬ 
dex directions. Note that, instead of using polar coordi¬ 
nates, almost the same result could be obtained by di¬ 
viding the horizontal width by the mean scattering angle 
directly from raw images. The more rigorous, polar ap¬ 
proach was adopted here with the advantage that the tilt 
angle can be measured directly. Fig. [6] shows the smooth 



Figure 5: Azimuthal scan of the polar width as measured by a 
Lorentzian fit (see e.g. Fig. HI). The scan outlines sharp varia¬ 
tions taking place within one or two degrees. The six labels (a) to 
(f) refer to the images in Fig. [2] 


3.1.3. ID scattering and polar profiles 
The scattering profile in the plane of incidence is compar¬ 
atively broad, and it is reasonably well fitted by a log¬ 
normal distribution commonly used in atom surface scat¬ 
tering EH. This profile is used here only to measure the 
position of the maximum. All polar plots P{a) exploited 
for the azimuthal scan correspond to a slice of -0.2° below 
the maximum of this scattering distribution. Two exam¬ 
ples of the polar angle distribution are displayed in Fig. HI 
and are empirically fitted by a Lorentzian profile to extract 
three parameters: the mean value (a), the width Wa, and 
the intensity. 



evolution of the mean polar angle {a)bottom measured by 
the fit on the bottom part of the scattering distribution 
that amounts to 0.2° below the maximum determined from 
the log-normal fit of the vertical profile. This dependence 
is adjusted with a pure sine function with a full amplitude 
of 4.3°. This corresponds to a tilt angle of r = 1.1° as 



Figure 4: Polar scattering profiles P(q;). The left panel (a) corre¬ 
sponds to the images displayed in Fig. He). while the right panel 
(b) is derived from Fig. Of), shown also in Fig. Qab 


Figure 6: Azimuthal scan of the center (o') bottom of polar profile 
(o) that is measured by the Lorentzian fit in Fig. [5]in the bottom part 
of the images. The diamonds (♦) plot the difference with the same 
quantity (a)top measured in the upper part of the same scattering 
images. 


sketched in Fig. [71 It should be mentioned that in the 
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natural angular frame of the detector (referred to as the 
Cartesian frame above) this 4.3° amplitude of the angular 
swing is only 0.05° since the typical value of Oeff repre¬ 
senting the radius of the circle in Fig. [3] corresponds to 
0.7°. As a naive attempt to capture the rotation of the 
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Figure 7: Sketch of the maximum effect of a tilt angle r between 
the surface (dark gray) and the rotation axis (light grey) on the 
specular reflection in the polar coordinates. Note that the main 
velocity component of the beam is not visible. 


ellipse visible in Fig. [2] we have also measured the mean 
polar angle {a)top^ a 0.2° wide slice on top of all images, 
and subtracted it from the earlier defined (a)bottom , he. 
the 0.2° slice below the maximum. Fig. [H] shows bipo¬ 
lar oscillation around the low index directions identified in 
Fig. [5l This nicely illustrates the rotation effect of the 
quasi-ellipse. 

Combining the azimuthal scan in Fig. [5] together with the 
initial GIFAD diffraction and FEED data is comparatively 
easy because each piece of information is complementary, 
however, this will not be discussed here. The reason is that 
we lack the redundancy of a complete azimuthal scan and, 
more importantly, the rapid degradation of the sample in 
a single day at a pressure below 10“^^ mbar. The lat¬ 
ter incentive demonstrates that the molecular organization 
is metastable reducing our confidence that the diffraction 
and triangulation measurements were actually observing 
the same structure. This doubt was also triggered by the 
fact that our observation does not directly match the di¬ 
rect STM imaging after a comparable preparation. This is 
probably due to the fact that the system seems very sen¬ 
sitive to the density of deposited molecules so that there 
could be a gradual transition from a densely packed com¬ 
paratively well-organized layer to a looser organization [l[ . 


4. Discussion 

The present experimental conditions were chosen to track 
diffraction features. This has led to meticulous recording 
of more than 100 successive images each exposed for 100 
seconds with additional manual control, stretching over al¬ 
most ten hours. One of the reasons for that was that the 
beam intensity was low, the other being due to the cor¬ 
rections imposed by the tilt angle. However, even without 
any tilt of the surface, i.e. without any need for man¬ 
ual control, this would have required three hours of beam 
time. The experiment duration can be drastically reduced 
by noting that the smaller structure observed in the im¬ 
ages is of the order of 0.08° full width at half maximum 


(EWHM) so that, if triangulation only had been targeted, 
the two 100jam diaphragms reducing the beam divergence 
below 0.02° could have been replaced by two 400jam di¬ 
aphragms giving an intensity 256 times larger. Indeed res¬ 
olution scales with the diameter of the diaphragm but the 
beam intensity scales with the fourth power of the diame¬ 
ter. A scan would take less than 10 min without any loss. 
It is possible to gain another order of magnitude to bring 
the acquisition time close to a minute by compromising 
the resolution or by using a more sophisticated image pro¬ 
cessing algorithm. Scanning times shorter than a minute 
are probably not relevant since fast mechanical action un¬ 
der vacuum is usually not recommended. These advances 
are still compatible with in-situ analysis given the ben¬ 
efit that molecular growth is usually more uniform on a 
rotating sample smearing possible beam inhomogeneity. 
Einally, the comparatively low angular resolution of the 
present study can be improved by a continuous recording 
while the image statistics can be recovered by averaging 
successive images. 

The width method detailed here is similar in essence to 
that reported in ^ where the peak intensity of the scat¬ 
tering profile was used. Both approaches can probably be 
further improved and generalized by a generic multipo¬ 
lar expansion of the recorded images to extract the main 
scattering parameters: mean values, width, peak inten¬ 
sity, and more refined ones like the k^-ky correlation terms 
quantifying the orientation of the scattering pattern. The 
other aspect illustrated here is that triangulation is less 
demanding than diffraction. This can be extremely useful 
to extend the range of GIEAD operation. Such a toler¬ 
ance to local disorder can be understood in terms of the 
surface coherence length Ls defined as the mean distance 
between defect. This distance should be compared with 
the mean length Lt of the atom trajectory on the sur¬ 
face estimated from classical trajectories or simple models 


2l|. Elastic diffraction can only occur if Ls ^ Lt and 


this has not been observed so far for molecular layers. If 
Ls ^ Lt inelastic diffraction can be observed providing 
accurate lattice parameter, but its quantitative interpre¬ 
tation in terms of surface electronic density might be dif¬ 
ficult. If Ls ^ Lt then no diffraction should be observed. 


however the momentum transferred to the projectile atom 
along its trajectory is still very much influenced by lo¬ 
cal molecular alignment as probed in triangulation. In 
other words triangulation is likely to identify the local or¬ 
der over a much broader range of surface coherence length 
(i.e. global order). Meanwhile inelastic diffraction and 
possibly elastic diffraction will progressively show up only 
as the quality of organic layer improves. Both IBT and 
GIEAD have demonstrated a sensitivity to the topmost 
layer and are able to count the exact number of succes¬ 
sive layers (HI, |23| during growth. Both approaches are 
fully compatible and obviously needed to help monitoring 
in-situ the growth of very well organized molecular lay¬ 
ers. The very encouraging new aspect demonstrated here 
is that azimuthal scans are well-contrasted even with the 
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low perpendicular energies needed for GIFAD. These low 
values, ^0.1 eV should not induce any damage to the 
fragile organic film. At first sight, it seems rather illusory 
to record an azimuthal dependence at one degree of inci¬ 
dence of a surface itself tilted by more than a degree! The 
straightforward image processing detailed here shows that 
an azimuthal scan with significant intensity variations can 
be retrieved. 

5. Conclusion 

Using atomic beams of very low effective interaction en¬ 
ergy with the surface in a regime where the damage to 
fragile molecular layers is expected to vanish constitutes 
a powerful method to probe the self-assembly of adsorbed 
monolayers. A simple processing of the scattering profile 
generates highly contrasted azimuthal scans revealing the 
directions of molecular alignment. In this respect, imple¬ 
menting triangulation with a position-sensitive detector is 
certainly a good choice paving the way to atomic diffrac¬ 
tion if a molecular layer with long-range order is formed. 
Additionally, triangulation becomes compatible with GI¬ 
FAD offering fewer restrictions on the level of organization 
so that molecular layer growth can be monitored over a 
very broad range of relative disorder. 
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